*Create Roc Plots*
* Use "FPA In Sample repdata.dta"for this figure*
*Hijacking Plot*
logit      hijack_dum ucdp_lag  gdpln_lag sfi_lag  wto_lntradelag polity2lag cinc_lnlag coast_ln  pop_ln t1hij t2hij t3hij if year<2007, cluster(ccode)
lroc, title("Hijackings") ytitle("True Positive Rate") xtitle("False Positive Rate") msize(vtiny) name(Hijackings)
*Incident Plot*
logit      pirincidents_dum  gdpln_lag ucdp_lag sfi_lag  wto_lntradelag polity2lag cinc_lnlag coast_ln  pop_ln  t1inc t2inc t3inc if year<2007, cluster(ccode)
lroc, title("Incidents") ytitle("True Positive Rate") xtitle("False Positive Rate") msize(vtiny) name(Incidents)
*Combine the three plots*
gr combine Incidents  Hijackings
